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1. Introduction 


Cooperative techniques are promising solutions for cellular wireless systems to improve 
system fairness, extend the coverage and increase the capacity. Antenna array schemes, also 
referred as MIMO systems, exploit the benefits from the spatial diversity to enhance the link 
reliability and achieve high throughput (Foschini & Gans, 1998). On the other hand, 
orthogonal frequency division multiplexing (OFDM) is a simple technique to mitigate the 
effects of inter-symbol interference in frequency selective channels (Laroia et al., 2004). The 
integration of multiple antenna elements is in some situation unpractical especially in the 
mobile terminals because of the size constraints, and the reduced spacing does not 
guarantee decorrelation between the channels. An effective way to overcome these 
limitations is generate a virtual antenna-array (VAA) in a multi-user and single antenna 
devices environment, this is referred as cooperative diversity. The use of dedicated 
terminals with relaying capabilities has been emerging as a promising key to expanded 
coverage, system wide power savings and better immunity against signal fading (Liu, K. et 
al., 2009). 

A large number of cooperative techniques have been reported in the literature the potential 
of cooperation in scenarios with single antennas. In what concerns channel estimation, some 
works have discussed how the channel estimator designed to point-to-point systems 
impacts on the performance of the relay-assisted (RA) systems and many cooperative 
schemes consider that perfect channel state information (CSI) is available (Muhaidat & 
Uysal, 2008), (Moco et al., 2009), (Teodoro et al., 2009), (Fouillot et al., 2010). Nevertheless, to 
exploit the full potential of cooperative communication accurate estimates for the different 
links are required. Although some work has evaluated the impact of the imperfect channel 
estimation in cooperative schemes (Chen et al., 2009), (Fouillot et al., 2010), (Gedik & Uysal, 
2009), (Hadizadeh & Muhaidat, 2010), (Han et al., 2009), (Ikki et al. 2010), (Muhaidat et al., 
2009), new techniques have been derived to address the specificities of such systems. 
Channel estimation for cooperative communication depends on the employed relaying 
protocol, e.g., decode and forward (DF) (Laneman et al., 2004) when the relay has the 
capability to regenerate and re-encode the whole frame; amplify and forward (AF) 
(Laneman et al., 2004) where only amplification takes place; and what we term equalize and 
forward (EF) (Moco et al., 2010), (Teodoro et al., 2009), where more sophisticated filtering 
operations are used. 
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In the case of DF, the effects of the BR (base station-relay node) channel are reflected in the 
error rate of the decided frame and therefore the samples received at the destination only 
depend on the RU (relay node-user terminal) channel. In this protocol the relaying node 
are able to perform all the receiver’s processes including channel estimation and the point- 
to-point estimators can be adopted in these cooperative systems. However the situation is 
different with AF and Equalize-and-Forward (EF) which are protocols less complex than the 
DF. In the former case (AF), B>RU (base station-relay node-user terminal) channel is the 
cascaded of the BÈR and RU channels, which has a larger delay spread than the 
individual channels and additional noise introduced at the relay, this model has been 
addressed in (Liu M. et al., 2009), (Ma et al., 2009), (Neves et al., 2009), (Wu & Patzold, 2009), 
(Zhang et al., 2009), (Zhou et al., 2009). 

Channel estimation process is an issue that impacts in the overall system complexity reason 
why it is desirable use a low complex and optimal estimator as well. This tradeoff has been 
achieved in (Ribeiro & Gameiro, 2008) where the MMSE in time domain (TD-MMSE) can 
decrease the estimator complexity comparatively to the frequency domain implementation. 
In (Neves et al., 2009) it is showed that under some considerations the TD-MMSE can 
provide the cascaded channel estimate in a cooperative system. Also regarding the receiver 
complexity (Wu & Patzold, 2009) proposed a criterion for the choice of the Wiener filter 
length, pilot spacing and power. (Zhang et al., 2009) proposed a permutation pilot matrix to 
eliminate inter-relay signals interference and such approach allows the use of the least 
square estimator in the presence of frequency off-sets. Based on the non-Gaussian dual-hop 
relay link nature (Zhou et al., 2009) proposed a first-order autoregressive channel model and 
derived an estimator based on Kalman filter. In (Liu, M. et al., 2009) the authors propose an 
estimator scheme to disintegrate the compound channel which implies insertion of pilots at 
the relay, in the same way (Ma et al., 2009) developed an approach based on a known pilot 
amplifying matrix sequence to improve the compound channel estimate taking into account 
the interim channels estimate. To separately estimate BÈR and RDU channels (Sheu & 
Sheen. 2010) proposed an iterative channel estimator based on the expectation 
maximization. Regarding that the B>R and RU links are independent and point-to-point 
links (Xing et al. 2010) investigated a transceiver scheme that jointly design the relay 
forward matrix and the destination equalizer which minimize the MSE. Concerning the two- 
way relay (Wang et al. 2010) proposed an estimator based on new training strategy to jointly 
estimate the channels and frequency offset. For MIMO relay channels (Pang et al. 2010) 
derived the linear mean square error estimator and optimal training sequences to minimize 
the MSE. However to the best of our knowledge channel estimation for EF protocol that use 
Alamouti coding from the base station (BS) to relay node (RN), equalizes, amplifies the 
signals and then forward it to the UT has not been considered from the channel estimation 
point of view in the literature. Such a scenario is of practical importance in the downlink of 
cellular systems since the BS has less constraints than user terminals (or terminals acting as 
relays) in what concerns antenna integration, and therefore it is appealing to consider the 
use of multiple antennas at the BS improving through the diversity achieved the 
performance in the BR link. 

However due to the Alamouti coding-decoding operations, the channel B>RU is not just 
the cascade of the B>R and RDU channels, but a more complex channel. The channel 
estimator at the UT needs therefore to estimate this equivalent channel in order to perform 
the equalization. The derivation of proper channel estimator for this scenario is the objective 
of this chapter. We analyze the requirements in terms of channels and parameters estimation 


www.intechopen.com 


Channels and Parameters Acquisition in Cooperative OFDM Systems 117 


to obtain optimal equalization. We evaluate the sensitivity of required parameters in the 
performance of the system and devise scheme to make these parameters available at the 
destination. We consider a scenario with a multiple antenna BS employing the EF protocol, 
and propose a time domain pilot-based scheme (Neves, et al. 2010) to estimate the channel 
impulse response. The B>R channels are estimated at the RN and the information about the 
equivalent channel inserted in the pilot positions. At the user terminal (UT) the TD-MMSE 
estimator, estimates the equivalent channel from the source to destination, taking into 
account the Alamouti equalization performed at the RN. The estimator scheme we consider 
operates in time domain because of the reduced complexity when compared against its 
implementation in frequency domain, e.g. (Ribeiro & Gameiro, 2008). 

The remainder of this chapter is organized as follows. In Section 2, we present the scenario 
description, the relaying protocol used in this work and the corresponding block diagram of 
the proposed scheme. The mathematical description involving the transmission in our 
scheme is presented in Section 3. In Section 4, we present the channel estimation issues such 
as the estimator method used in this work and the channels and parameters estimates to be 
assess at the RN or UT. The results in terms of BER and MSE are presented in Section 5. 
Finally, the conclusion is pointed out in Section 6. 


2. System model 


2.1 Definition 

Throughout the text index n and k denote time and frequency domain variables, 
respectively. Complex conjugate and the Hermitian transposition are denoted by (.) and 
(-)", respectively. E{} and (*) correspondently denote the statistical expectation and the 
convolution operator. M (m,o*) refers to a complex Gaussian random variable with mean 
m and variance o°. diag(-) stands for diagonal matrix, || denotes absolute value and I, 


denotes the identity matrix of size Q . Regular small letters denote variables in frequency 
domain while boldface small and capital letters denote matrices and vectors, respectively in 
frequency domain as well. Variables, vectors or matrixes in time domain are denoted by 
(~). All estimates are denoted by (^). 


2.2 Channel model 

The OFDM symbol x=d+p where p corresponds to the pilots which are multiplexed 
with data d subcarriers. The element x, of the OFDM symbol vector is transmitted over a 
channel which the discrete impulse response is given by: 


hi = >, 85(n-t,), (1.1) 


where G is total number of paths, 2, and r, are the complex amplitude and delay of the 
gth path, respectively. Ø, is modelled as wide-sense stationary uncorrelated scattering 
(WSSUS) process. The gth path has a variance o? which is determined by the power delay 
profile and satisfies 2: o, =1. Although the channel is time-variant we assume it is 
constant during one OFDM symbol interval and its time dependence is not present in 


notation for simplicity. 
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2.3 Scenario description 

The studied scenario, depicted in Fig. 1, corresponds to the proposed RA schemes for 
downlink OFDM-based system. The BS and the RN are equipped with M and L antennas, 
respectively. The BS is a double antenna array and the UT is equipped with a single 
antenna. Throughout this chapter we analyze two RA schemes: the RN as a single antenna 
or an array terminal. These scenarios are referred as MxLx1 schemes. 
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a ru 
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Ô - - -. Relay Channel 
Direct Channel 
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Fig. 1. Proposed RA scenario 


The following channels per k subcarrier are involved in this scheme: 
e 7 =Mx1 MISO channel between the BS and UT (B>U): h ,m=1,2 


bum,(k) 


e Lx1 MISO channel between the RN and UT (R9U): h ,1=1,2 


rul,(k) 
e MxL channel between the BS and RN (BOR): h 
All the channels are assumed to exhibit Rayleigh fading, and since the RN and UT are 


mobile the Doppler’s effect is considered in all channels and the power transmitted by the 
BS is equally allocated between the two antennas. 


,m=1,2 andl=1,2 


brml,(k) 


2.4 The Equalize-and-Forward (EF) relaying protocol 

For the single antenna relay scenario, the amplify-and-forward protocol studied in (Moco et 
al., 2010) is equivalent to the RA FF protocol considered here. However, if the signal at the 
relay is collected by two antennas, doing just a simple amplify-and-forward it is not the best 
strategy. We need to perform some kind of equalization at the RN to combine the received 
signals before re-transmission. Since we assume the relay is half-duplex, the communication 
cycle for the aforementioned cooperative scheme requires two phases: 

Phase I: the BS broadcasts its own data to the UT and RN, which does not transmit data 
during this stage. 

Phase II: while the BS is idle, the RN retransmits to the UT the equalized signal which was 
received from the BS in phase I. The UT terminal receives the signal from the RN and after 
reception is complete, combines it with the signal received in phase I from the BS, and 
provides estimates of the information symbols. 
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2.5 The cooperative system 

Fig. 2 shows the corresponding block diagram of the scenario depicted in Fig. 1, with 
indication of the signals at the different points. The superscripts (1) and (2) denote the first 
and the second phase of the EF protocol, respectively. In the different variables used, the 
subscripts u, r and b mean that these variables are related to the UT, RN and BS, 
respectively. 


1) 2(1) 
7 AN Oo ) 
>. Data Hard 4 
Soft-Decision Combiner | | Decision |, : 
pe i (2) o" 
$ Su(k) A 
k i Estimator i 
id == iat 


er eee 


Fig. 2. The corresponding block diagram of the MxLx1 RA scenario 


Let d=(d, dy +++ dy, =i). be the symbol sequence to be transmitted where N, is the 
number of data symbols, then for k even the SFBC (Teodoro el al., 2009) mapping rule is 
defined in Table 1. The symbols d, are assumed to have unit average energy, i.e. 
[do 
energy transmitted by the two antennas per subcarrier is normalized to 1. 


=1, Vk, and therefore the factor 1/ V2 used in the mapping, is to ensure that the total 


Subcarrier Antenna 1 Antenna 2 


Table 1. Two transmit antenna SFBC mapping 


The pilot symbols are multiplexed with data and the BS broadcasts the information ras 
(data and pilot) to the RN and UT. This processing corresponds to the phase I of the EF relay 
protocol. At the UT, the direct channels are estimated and the data are SFBC de-mapped and 
equalized. These two operations are referred as soft-decision which the result is the soft- 
decision variable, in this case, Ce At the RN, pilots and data are separated; based on 
pilots, the channels BƏR are estimated and the soft-decision is performed. The result is the 
soft-decision variable se Then, the new pilot symbols are multiplexed in Su and the 
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information a is transmitted / forwarded to the UT via R>U channel. This second 
transmission corresponds to the phase II of the EF protocol. At the final destination, the 
required channel is estimated and the soft-decision is performed in order to obtain the soft- 
decision variable ca . After the phase II the UT has the soft-variable provided by both the 
BS and RN. These variables are combined and hard-decoded. 


3. Mathematical description of the proposed cooperative scheme 


The mathematical description for transmit and receive processing is described in this 
section. As this work is focused on channel estimation, this scheme is designed in order to 
be capable to provide all the channels and parameters that the equalization requires in both 
phases of the relaying protocol. 


3.1 Phase | 
During the first phase the information is broadcasted by the BS. The frequency domain (FD) 
signals received at the UT in data-subcarriers k and k+1 are given by 


(cde B . (1) 
Yun = No Cea Dara] + 1y (x) 
j (3.1) 
1 1 * 1 
i ai z luo + E T E 
where i y 18 the additive white Gaussian noise with zero mean unit variance ol) and for 
m=1,2, hum represent the channels between the BS and the UT terminal. 
The FD sionals received at the RN in data-subcarriers k and k+1 are expressed by: 
Yi = zl br11,( nets) i Th) 
„l = 1, 2 (3.2) 


Ya eae 2 te, ww) + E an) + Wy (e412) 


where h ati the channels between the antenna m of the BS and antenna / of the 


bril,(k) 
RN terminal and nl *) is the additive white Gaussian noise with zero mean unit variance 


ase 
Sind the data are SFBC mapped at the BS the SFBC de-mapping at the terminals RN and UT 
also includes the MRC (maximum ration combining) equalization which coefficients are 
functions dependent on the channels estimates. It is widely known that in the OFDM 
systems the subcarrier separation is significantly lower than the coherence bandwidth of the 
channel. Accordingly, the fading in two adjacent subcarriers can be considered flat and 
. Thus, in phase I 


without loss of generality we can assume for generic channel hiy = hgs) 


the soft-decision variables at the UT follow the expression. 


1 *(1 
k) = Sout (k otie 


() 
(3.3) 


u,(k) 
(1 “(1 * 1 
Sa (ket) = Sae T Pity a 
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where the equalization coefficients for m=1,2 are given by gpm (4) = 


»/v2 V2 o.. After 


some mathematical manipulation, these soft-decision variables may be expressed as: 


Pam, ( 


hy h 

(1 bul,(k) (1) bu2,(k) (1) 
sn, =r d+ n+ n 

w(k) bu,(k) (k) / u,(k) u,(k) 

E ie (3.4) 
h, h 

Q bul(k+1) (1)  _  bu2,(k+1) (1) 

Su,(k+1) = Duet) est) + J2 nN, (k+1) 2 Na (k+1) 


2 


bum 


1% 
where Ta a h 


m=1 
The soft-decision variables should be kept in a buffer, waiting for the information to be 
provided by the RN in the second phase of the protocol. The mathematical formulation of 
the next phase varies according to the number of antennas at the RN, i.e. L, and these cases 
are separately presented in the next sub-sections. 


3.2 Phase Il 
3.2.1 RN equipped with a single antenna 
The FD soft-decision variables at the RN in data-subcarriers k and k+1 are expressed by: 


ys Tye + Baten 
1) 


4! 
M i (3.5) 


S,,(k+1) 7 n, (k) RTA k+1) 


where the equalization coefficient are given by Sy. (x) = Fat / V2 ) ,for m=1,2 and l=1 


After some mathematical manipulation, these soft-decision variables are given by: 


(3.6) 


=> 


&®@ -=r d { hya) (1) “br 21, (k+1) o 


r,(k+1) br,(k+1)”(k+1) J2 r,(k+1) OAR 


Tair 
where Tor) = 5 hypen 


m=1 
In order to transmit a unit power signal the RN normalizes the expression in (3.6) by 
considering the normalization factor @,) which is given by: 


1 
ay =- (8.7) 


2(1 
Fao + Tg. ) 

During the phase II, the normalized soft-decision variable is sent via the second hop of the 

relay channel RU. The FD signal received at the UT per k subcarrier is expressed 

according to: 
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(2) _ 1) 2) 
Yau (k) = XoSn iu) Tua (3.8) 


where ny ly is the additive white Gaussian noise which is zero mean and has unit variance 
2(2 
ož ) 


The signal in (3.8) is equalized using the coefficients Seat) = (a0 m / o?) which after 
some mathematical manipulation is given by: 


2 


2) _ 272 rul,(k) rul(k) _ (2) 
8 =a, do +æ 0 n 
(k) (k) br,(k) o? (k) (K) (k) o? (k) (3.9) 
z 2 1 A 1) r 
+a, ie k rul,(k br11,(k nl k + fizi k n, k 
(k) br,(k) W fp z| (K) (hk) (k) D) 


The equalization coefficient g is a function dependent on the channel estimate f 


rul,(k) rul,(k) 
and the variance of the total noise. Moreover, the statistics of the total noise is conditioned to 


the channel realization h ,)- Therefore the variance of the total noise can be computed as 


brml,( 
conditioned to these channel realizations or averaged over all the channel realizations. We 


denote by On 


w) the noise variance conditioned to the specific channel realization per k 
subcarrier and by o; the unconditioned noise variance. The noise variance of the total noise 


conditioned to channel realizations is found to be: 


2 2(1 
E inaa = T br (4) Pran) Pea (3.10) 


3.2.2 RN equipped with a double antenna array 
When the array is equipped with two antennas the soft-decision variables are expressed by: 


2 
1 1 1)" 

ah = Faa + S a) 

p ,l=1,2 (3.11) 

1 1)* 1 

AOR X VA Eont ety + MATEN 


{}=1 


where the equalization coefficients are expressed by 8n) = (lop / V2 ) ‘ 


The RN transmits to the UT a unit power signal following the normalization factor in (3.7). 


A 2 
Pita] . The FD signals 


. . 1 2. o2 
However, in this scenario Tœ) is expressed by I’, a ap 


m=1 l=1 


received at the UT are given by: 


eee O _ = a0 (2) 
Yuk) = Ca hoita a] HALIK) 
A , (3.12) 
(2 _ ees)" 1) (2) 
Wea = —=(h a, rae +h OTE ) +n 


rul,(k+1)~(k+1)°r,(k+1) u,(k+1) 
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where Taig represent the channels between the RN and the UT terminal. The soft-decision 


variables found at the UT in phase II are expressed by: 
(2) Zo ye *(2) 
a(k) 7 Bai F Sru2Yu,(k+1) 
2 *(2 7 2 
ay = -8Y + Bea ta 


S 
' (3.13) 


where the equalization coefficients are given by 8,4%) = (aol cedar i V20?) , for l=1,2. 
In this scenario the soft-decision variables also depend on the variance of the total noise of 
which conditioned to the channel specific realization can be calculated from (3.12) and is 
expressed by, 


5 a aaa ee (3.14) 


Ip p 
where lw) = m2 h 


rul,(k) 


The UT combines the signals received from the RN and the BS. By performing this 
combining the diversity of the relay path is exploited. This processing is conducted by 
taking into account Ea + a with ensure and MRC combining. The result corresponds to 
the variable to be hard-decoded. 


4. Channel estimation 


4.1 Time Domain Minimum Mean Square Error (TD-MMSE) estimator 

The TD-MMSE (Ribeiro & Gameiro, 2008) corresponds to the version of the MMSE estimator 
which was originally implemented in FD. This estimator comprises the least square (LS) 
estimation and the MMSE filtering, both processed in time domain (TD). The TD-MMSE is a 
pilot-aided estimator, i.e. the channel estimation is not performed blindly. It is based on pilot 
symbols which are transmitted by the source and are known at destination. 

The pilot subcarriers convey these symbols that are multiplexed with data subcarriers 
according to a pattern, Fig. 3, where N f and N, correspond to distance between two 
consecutives pilots in frequency and in time, respectively. N is the number of OFDM 
symbols and N, is the number of subcarriers. The pattern presented in Fig. 3 is adopted 
during the transmission stages of the envisioned cooperative scheme. 

It is usual the pilot symbols assume a unitary value and be constant during an OFDM 
symbol transmission. Thus, at k subcarrier the element p of the vector p may be expressed 
by a pulse train equispaced by N, with unitary amplitude. The corresponding expression 
in TD is also given by a pulse train with elements in the instants (n-mN, iN), for 
me {0 iN; - 1} , according to the following expression. 


= ae ae 
Pr) E Paes Pin) E N. Ol enini] (4.1) 


The transmitted signal is made-up of data and pilot components. Consequently, at the 
receiver side the component of the received signal in TD is given by the expression in (4.2). 
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In order to perform the LS estimation in TD, i.e. TD-LS, the received signal is convolved 
with the TD pilot symbols ñ, . This convolution corresponds to multiply by 1 the 
subcarriers at frequencies N, as by design these are the positions reserved to the pilots thus 
the data component in the received signal vanishes. The resulting CIR estimate h,, is made- 
up of N, replicas of the CIR separated by N./N; . 


Aa 


x Nya. Ny-1 N 
i 2 = 0/1 Ne 1. 4. 
his 2 hamn) F D Mn, JN) 7 0, recy N, ( 3) 


CIR Noise 


Besides to estimate the CIR the TD-MMSE in (Ribeiro & Gameiro, 2008) can estimate the 
noise variance as well. It corresponds to an essential requirement when the UT has no 
knowledge of this parameter. Since we know that the CIR energy is limited to the number 
of taps, or the set of the taps {G} , the noise variance estimate ôf can be calculated by take 
into account the samples out of the number of taps, i.e. n¢{G} and by averaging the 
number of OFDM symbols N . Thus ĉ? is given by: 


2 


hn) . 


a2 
O, 


N,/N 
NN )-G]N 22 


(4.4) 


The MMSE filter can improve the LS estimates by reducing its noise variance. The TD- 
MMSE filter corresponds to a diagonal matrix with non-zero elements according to the 
number of taps, i.e. G , thus it can be implemented simultaneously with the TD-LS estimator 
and this operation simplifies the estimator implementation. The MMSE filter implemented 
by the (N ; /N f ) x (N -/N s) matrix and for a generic channel h itis expressed by: 
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Wumse,i = RaRa (4.5) 
where R., is the (N -/N 7 )x(N, /N i) filter input correlation matrix, E! hh” t, which is given 
by R,, +o ly, a and R- is the (N, /N Jr N. /N a] filter input-output cross-correlation 
matrix, Efhh |, which is given by R; = diag ( Oi Op; Obar Op 0] 

4.2 Channels and parameters estimates 

According to the scenario previously presented, there are channels which correspond to 
point-to-point links: Apan (x) and Mym- Therefore, these channels can be estimated by using 
conventional estimators. However, for the RN>UT links, and since the EF protocol is used, 
it is necessary to estimate a version of Hawa which depends on A 4) and Pay, the 
equivalent channel h. = @ 0 (yux) Note that the UT has no knowledge of a, and Tù. 
These factors are dependent on h.g which the UT has no knowledge as well. However, 
the channels hmg are estimated at the RN, and based on that, aI) is calculated. 
Therefore, we propose to transmit the factor a,I,, at the pilot subcarriers as pilots. 
Consequently, the new pilots are no longer constant and that may compromise the 
conventional TD-MMSE performance, since this estimator was designed in time domain 
assuming the pilots are unitary with constant values at the destination. Although our 
approach enables the destination had knowledge of the non-constant pilot (alw) , the 
result of the convolution between the received signal and these pilots, would result in the 
overlapped replicas of CIR, according to the Fig. 4. 
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Frequency Domain Time Domain 
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constant and undesirable 
samples between them 
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Equispaced pilots with 


Data samples overlap with 
non-constant values 
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Time Domain 
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Convolution by the pulse 
are spread 


train in TD 


Fig. 4. Pilots with non-constant values result in the overlapped replicas of the CIR 


However, the @,, expressions depend on the noise variance of ®© and the product liy 
tends to one for a high SNR value, according to (4.6). The same equation also suggests that 
the factor œl x) varies exponentially according to the SNR, as depicted in Fig. 5 for L=1. 


1 1 
OAN BE T =| — 
(K) (k) 2 br,(k) 2 
ae + Dy.) a Tir + 0 


hayek (4.6) 
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Fig. 5. QI vs. SNR 


Other behaviour of the factor a, ry can be demonstrated in terms of SNR and subcarriers, 
as presented in Fig. 6. These plots show that in the first case, i.e. SNR =20dB, the a,F() 
factor presents the amplitude close to 1 with some negligible fluctuation. However, in the 
second case, SNR = 2dB , the result is likely different to the previous one: the AT ay factor 
presents an amplitude also close to 1 but the fluctuation is not negligible. 
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Fig. 6. Qali per subcarriers 


The results in Fig. 6 emphasize that transmit the factor æl} in the pilot subcarriers may 


(k) 
degrade the estimator performance and the causes are: 
1. Pilots with some fluctuation in amplitude: 
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e As the amplitude of the pilots at the destination are constant and equal to one, the 
result of the estimation is a spread of the replicas of the CIR. 

2. Decreasing the amplitude of the pilots: 

e The SNR of the pilots is decreased as well. 
3. The MMSE filter depends on the statistics of the channel BPR 
Despite we are considering the TD-MMSE estimator in our analysis, the causes presented 
previously degrade the performance of any other estimator scheme as well. In order to 
quantify how these effects can degrade the estimator performance we evaluated the impact 
of both of them, separately, in a SISO system, since the B>R and R>U channels correspond 
to point-to-point links. 
First, we evaluate the case when the pilots have some fluctuation in amplitude. We consider 
that the pilots (originally with unit i ea had their amplitude disturbed by a noise 
with zero mean and variance equal to oĉ, =Ef1- Aal) | i. oĉ- quantifies how far the 


factor ælg would be from the pilots with unitary a. We can express oĉ, as: 


o? =1+ el( oul a) ï) -Æla Ty}: (4.7) 


Therefore, the pilots are no longer constant and unitary. They have some fluctuation in 
amplitude which depend on aI, and they are eau to Pa = 1+z, where z has a 
normal distribution with zero mean and power o.,. The performance of a SISO system 
which the pilots correspond to p, 2 Is shown in Fig, 7 (dash line). For reference, we also 
include the SISO performance for unit pilots, p,. Since we are focus on the degradation of 
the estimator performance, the results are presented in terms of the normalized mean square 
error (MSE) and E,/N,, where E, corresponds to the energy per bit received at UT and N, 
is the power spectrum density of the total noise which affects the information conveying 
signals. The normalized MSE for a generic channel h is given by: 


s] T 


For low values of SNR, [0 — 4] , the major difference in performance between the two results 
is approximately 0.5dB, which is not a noticeable degradation and the estimator 
performance is not compromised. 

The second effect to be evaluated is the decreasing of the amplitude of the transmitted pilot. 
In order to evaluate this effect we also consider a SISO system, for which the transmitted 
pilots assume different constant values, and we consider different values of the factor 
Al as pilot. The normalized MSE is given by (4.9). The result is shown in Fig. 8 and for 
reference we include the SISO performance for unit pilots, p,,as well. 


= elisi] (4.9) 
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Fig. 8. Channel estimation MSE performance 


The results show a constant shift in the MSE when the amplitude of the pilots is not unitary. 
The shift presents in all results is not a real degradation. It is caused by the normalization 
present in the MSE in (4.9). In fact, assuming a MSE without normalization the results are all 
the same. 


Transmit the factor aI) as pilot does not bring any noticeable degradation in the TD- 


MMSE performance eae to transmitting unitary pilots. The major degradation occurs 
only when the pilots have some fluctuation in amplitude, as shown in Fig. 7, and solely for 
low SNR values. 

The conventional MMSE filter in (4.5) is implemented to improve a channel estimate h 


when the required channel corresponds to h. However, according to our cooperative 
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scheme, we need estimate an equivalent channel h,, = a, Tal eut,(&) Since the factor al does 


k 


not depend on Å, , the MMSE filter input correlation matrix for the channel hi, , referred 


ru / eq / 


2 . é è 
~R,- ,- +I while the filter input-output 
Negheg ba Cae arar) hruħru Y: No/N¢ P P 


R. . , is expressed by efh ha) = Risa R 


cross-correlation, referred as R. . , is given by Efi h” | “Rea R- - , both R, . and 


heqħeq a arar| hruħru heqheq 
R. ; are (N -/N s)* (N -/N ‘) matrices. Thus the MMSE filter, when he is required, may 
teqeq 


be express as: 


-1 
_ 2 
Wass, feq RTE iie PRR is + oa . 


(4.10) 
As we shown previously in (4.6) the factor œI tends to one for high values of SNR and 
examining (4.10), which depends on ar , it is clear that (4.10) tends to (4.5) for high values 
of SNR as well. In order to show that several simulation were performed by taking into 


account R. i and the noise variance ol), According to Fig. 9 the results show that the 
maximum value in the R, . matrix is close to -40dB for high values of the noise variance. 
teq'teq 
ii | | | eae 
BG ljhs a renee seer ee yee | 
| | H + | 
| eeeenini ee ee a 
L | + | | 
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Fig. 9. Maximum value in the correlation matrix vs. noise variance 


According to the results in Fig. 7 and Fig. 8 in terms of MSE, transmitting the factor al 
brings, in the worst case, 0.5dB of degradation reason why there is no need to increase the 
system complexity by implementing the filter in (4.10). We have shown that our cooperative 
scheme allows the use of the TD-MMSE estimator without compromising its estimate. 
According to (4.6) the behaviour of the factor æl does not depend on number of taps of 


the channel, it depends only on the noise variance o2” . Therefore, this analysis is applied to 
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any other channel without loss of generality. Besides to estimate the equivalent channel it is 


2 
hag and ofl yl mqy for L=1 and L=2, 


ky br,(k)” ru, 


necessary estimate others factors ail 


respectively. These factors are required parameters in the variance of the total noise 


y previously presented in (3.10) and (3.14). 
) 


t (Mpa fk 
Although the UT does not have individual knowledge of hoem)” OL) and T (K) it has 
knowledge of the second moment of the expected value of the channels, i.e. for all channels 


Ela} =1. Thus, we propose the use of the noise variance unconditioned to the channel 


e . T . : . 2 2 
realization, o; , instead of its instantaneous value Oh F Therefore, of is referred as the 


rul,(k) 
expectation value of the variance of the total noise. Also we consider that the channels have 
identical statistics, i.e. ol) = ol) 2( 


and (4.12) for L=1 and L =2 , respectively. 


=o.) thus o? can be expressed numerically by (4.11) 


1 
Dens 2(2) 2(2) 
o = isso oho, (4.11) 
1 
Dns 2(2) 1 (2) 
o = TEE 20° +0. (4.12) 


If we consider the premise of the cooperative transmission, high SNR compared to the SNR 


2(2) 


of the direct link, we have o4“ <<1.5 for L=1 and consequently (4.11) may be express as 


oz SPO, 
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Fig. 10. System performance 
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To assess the validity of using the averaged noise variance instead of the conditioned one 
we plot in Fig. 10, the BER versus E,/N, performance assuming perfect channel estimation 
is available at the receiver but considering the cases where the noise variance used is the 
conditioned one and the averaged ones. The results refer to a channel as referred in Section 5 
but similar results were obtained with other models. 

The performance penalty by using the averaged noise variance is less than 0.8dB which is a 
tolerable penalty to pay in order to obtain the variance of the total noise regarding the low 
complexity implementation. Therefore we consider the use of o? in our schemes. 


5. Results 


5.1 Simulation parameters 

In order to evaluate the performance of the presented RA schemes we considered a typical 
scenario, based on LTE specifications (GPP TS, 2007). In the simulations we used the ITU 
pedestrian channel model B at speed v=10km/h. The transmitted OFDM symbol carried 
pilot and data with a pilot separation N, =4 and N, =1. 

We focus our analysis on the 2x1x1 and 2x2x1 scenarios and the simulations were 
performed assuming that the channels are uncorrelated, the receiver is perfectly 
synchronized and the insertion of a long enough cyclic prefix in the transmitter ensures that 
the orthogonality of the subcarriers in maintained after transmission. We use the TD-MMSE 
to estimate all the noise variances as well. 

The results are presented in terms of BER and MSE, both as function of E,/N,. The 
normalized MSE is defined according to (4.9). The MSE performance of the cooperative 
channel is evaluated by averaging the MSE’s of the direct and the relaying channel (Kim et 
al., 2007). Since the direct channel corresponds to a MISO its MSE is obtained also by 
averaging the MSE of the BU channels, both normalized. The MSE of the relaying channel 
corresponds to the MSE of the equivalent channel h = Alh g Which is calculated 


according to (5.1). Thus the resulting MSE, i.e. the MSE of the cooperative channel, is given 
by: 


MSE = 1 ( : (MSE, )+MSE,,_ } (5.1) 


5.2 Performance evaluation 

In order to validate the use of the proposed scheme, some channel estimation simulations 
were performed using the TD-MMSE estimator. Fig. 11 depicts the BER attained with 
perfect CSI and the TD-MMSE estimator when the RN was employing the proposed pilots. 
The difference of performance is minimal in most of the cases and in the 2x1x1 scheme 
which is in the worst case this difference is 0.5dB . 

Fig. 12 depicts the normalized MSE’s performance of the 2x11 scheme. These results show 
that the proposed pilot allocation, at the RN, according to Section I.B, allows the TD-MMSE 
satisfactory estimate the required channel. When comparing the channel estimator for the 
link with relay against the one of the direct link, there is some penalty which accounts for 
the additional noise added at the relay. The relative penalty decreases as E,/N, increases 
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and can be verified to converge to 2.2dB which is the factor of 5/3 that relates the total and 
individual noises in the asymptotic case of high SNR. According to Fig. 13, this penalty is 
smaller in the 2x21 scheme, since the factor aE Gy presents a flatter behavior. 


ToT TT 


E,/Npo of the direct link (dB) 


Fig. 11. System performance: RA 2x11 and RA 2x21 schemes 
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Fig. 12. Channel estimation MSE performance: RA 2x11 scheme 
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Fig. 13. Channel estimation MSE performance: RA 2x21 scheme 


6. Conclusion 


In this chapter we considered two problems of channel estimation in a scenario where 
spatial diversity provided by SFBC is complemented with the use of a half-duplex relay 
node employing the EF protocol. The channel estimation scheme was based on the TD- 
MMSE which led to a significant complexity reduction when compared to its frequency 
domain counterpart. We proposed a scheme where the estimates of the B>R link are 
inserted in the pilot positions in the ROU transmission. For the estimation of the equivalent 
channel, i.e. BRU, at the destination we analyzed several simplifying options enabling 
the operation of channel estimation namely the use of averaged statistics for the overall 
noise and the impact of the fluctuations in the amplitude of the equivalent channel. In the 
RA 2x1x1 scheme is shown that in the asymptotic case of high SNR, and equal noise 
statistics at the relay and destination the penalty in the estimation equivalent channel is 
2.2dB relatively to the case of a direct link using the same pilot density. This difference in 
performance is smaller in the RA 2x2x1 scheme since the equivalent channel presents a 
flatter behaviour. The resulting estimation was assessed in terms of the BER of the overall 
link through simulation with channel representative of a real scenario and the results have 
shown its effectiveness despite a moderate complexity. 
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